Color-based feature identification

ABSTRACT

Disclosed are embodiments of systems and methods for identifying features using color information in an image. The image may be formed from one or more display images comprising color information and features or feature components. Because color information may be used to identify features, more than one feature or feature component may be displayed in a display image. Because a plurality of features may be identified in a calibration image, an image system, such as a projector-camera system, can reduce the number of display images needed to calibrate the system.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation-in-part of and claims the priority benefit under 35 U.S.C. §120 of U.S. patent application Ser. No. 11/623,016, filed on Jan. 12, 2007 (pending), which claims the benefit of U.S. patent application Ser. No. 11/613,129, filed Dec. 19, 2006 (pending), which claims the benefit under 35 U.S.C. §119(e) of U.S. Provisional Patent Application No. 60/762,178, filed Jan. 24, 2006, each of which is incorporated herein by reference in its entirety.

BACKGROUND

1. Field of Invention

The present invention is related to methods of calibrating projectors in projector-camera systems, and more specifically related to methods of simplifying application of dual photography.

2. Description of Related Art

When projectors and cameras are combined, hybrid devices and systems that are capable of both projecting and capturing light are born. This emerging class of imaging devices and systems are known in the research community as projector-camera systems. Typically, the images captured by one or more cameras, is used to estimate attributes about the display environment, such as the geometric shape of the projection surfaces. The projectors in the system then adapt the images projected so as to improve the resulting imagery.

A key problem that builders of projector-camera systems and devices need to solve is the determination of the internal imaging parameters of each device (the intrinsic parameters) and the determination of the geometric relationship between all projectors and cameras in the system (the extrinsic parameters). This problem is commonly referred to as that of calibrating the system.

In the computer vision community, there is a large body of work for calibrating imaging systems with one or more cameras. A commonly used method in the computer vision community for calibrating cameras is described in article, “A flexible new technique for camera calibration”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 22 (11):1330-1334, 2000, by Zhengyou Zhang, which is herein incorporated in its entirety by reference. In this method, multiple images of a flat object marked with a number of known feature points (typically forming a grid) are captured by a camera, with the flat object posed at a variety of angles relative to the camera. The image location of each feature point is extracted, and since the relative location of each feature point is known, the collection of feature point locations can then be used to calibrate the camera. When two or more cameras are present in the system, the intrinsic parameters as well as the geometric relationship between all cameras can be estimated by having all cameras capture an image of the flat object at each pose angle.

Since projectors and cameras are very similar in terms of imaging geometry, it might seem reasonable that techniques for calibrating cameras may be applicable to the calibration of projectors, and projector-camera systems. However, since all camera calibration techniques require that the camera (i.e. the imaging device being calibrated) capture a number of images, it would appear that camera calibration techniques cannot readily be applied to projectors, which cannot capture images.

Therefore, in traditional systems the cameras are calibrated first, and then the calibrated cameras are used to calibrate the projector. In these systems, a “bootstrapping” procedure, in which a pair of calibrated cameras is used to form a stereo pair, is used. As it is known, a stereo pair of cameras can be used to estimate depth, i.e. a perspective view, of feature points visible to the stereo pair. To calibrate a projector, the projector is first made to project feature points onto a display environment (i.e. a projection surface), and then by using the stereo pair of calibrated cameras to determine the perspective depth location of the projected points, the projector can be calibrated accordingly to compensate for surface irregularities in the display environment.

While this bootstrapping technique is a tested-and-proven calibration method for projector-camera systems, it is not applicable to the calibration of self-contained projector-camera devices, since it requires the use of external stereo camera pairs.

A technique called dual photography was proposed by Sen et al. in article, “Dual Photography”, Proceedings ACM SIGGRRAPH, 2005, which is herein incorporated by reference in its entirety. Dual photography makes use of Helmholtz reciprocity to use images captured with real cameras to synthesize pseudo images (i.e. dual images) that simulate images “as seen” (or effectively “captured”) by projectors. That is, the pseudo image simulates a captured image as “viewed” by the projector, and thus represents what a projector-captured image would be if a projector could capture images. This approach would permit a projector to be treated as a pseudo camera, and thus might eliminate some of the difficulties associated with calibrating projectors.

Helmholtz reciprocity is based on the idea that the flow of light can be effectively reversed without altering its transport properties. Helmholtz reciprocity has been used in many computer graphics applications to reduce computational complexity. In computer graphics literature, this reciprocity is typically summarized by an equation describing the symmetry of the radiance transfer between incoming and outgoing directions, ω_(i) and ω_(o): ƒr(ω_(i)→ω_(o))=ƒr(ω_(o)→ω_(i)), where ƒr represents the bidirectional reflectance distribution function, BRDF, of a surface.

Thus, dual photography ideally takes advantage of the dual nature (i.e. duality relationship) of a projected image and a captured image to simulate one from the other. As is described in more detail below, dual photography (and more precisely Helmholtz reciprocity) requires the capturing of the light transport between a camera and a projector.

When dealing with a digital camera and a digital projector, dual photography requires capturing each light transport coefficient between every camera pixel and every projector pixel, at the resolution of both devices. Since a digital projector and a digital camera can both have millions of pixels each, the acquisition, storage, and manipulation of multitudes of light transport coefficients can place real practical limitations on its use. Thus, although in theory dual photography would appear to offer great benefits, in practice, dual photography is severely limited by its physical and impractical requirements of needing extremely large amounts of computer memory (both archive disk-type memory and active silicon memory) and needing extensive and fast computational processing.

A clearer understanding of dual photography may be obtained with reference to FIGS. 1A and 1B. In FIG. 1A, a “primal configuration” (i.e. a configuration of real, physical devices prior to any duality transformations) includes a real digital projector 11, a real projected image 13, and a real digital camera 15. Light is emitted from real projector 11 and captured by real camera 15. A coefficient relating each projected beam of light (from real projector 11) to a correspondingly captured beam of light (at real camera 15) is called the light transport coefficient. Using the light transport coefficient, it is possible to determine the characteristics of a projected beam of light from a captured beam of light.

In the present example, real projector 11 is a digital projector having an array of projector pixels 17 symbolically shown in a dotted box and comprised of s rows and r columns of projector pixels e. The size of projector pixel array 17 depends on the resolution of real projector 11. For example, an VGA resolution may consists of 640 by 480 pixels (i.e. 307,200 projector pixels e), an SVGA resolution may have 800 by 600 pixels (i.e. 480,000 projector pixels e), an XVG resolution may have 1024 by 768 pixels (i.e. 786,732 projector pixels e), an SXVG resolution may have 1280 by 1024 (i.e. 1,310,720 projector pixels e), and so on, with greater resolution projectors requiring a greater number of projector pixels e.

Similarly, real camera 15 is a digital camera having an array of light receptor pixels 19 symbolic shown in a dotted box and comprised of u rows and υ columns of receptor pixels g. The size of receptor pixel array 19 again depends on the resolution of real camera 15. However, it is common for real camera 15 to have a resolution of 4 MegaPixels (i.e. 4,194,304 receptor pixels g), or greater.

Since real projector 11 is a digital projector having an array of individual light projection pixels e and real camera 15 is a digital camera having an array of individual light receptor pixels g, a matrix T is used to describe the group of light transport coefficients relating each projector pixel e in real projector 11 to each receptor pixel g of real camera 15 (i.e. element T_(ge) in matrix T would be the transport coefficient from an individual, real projector pixel e to an individual, real camera receptor pixel g). Therefore, a real captured image C′ is related to a projected image P′ as C′=TP′.

The duality transformation, i.e. dual configuration, of the system of FIG. 1A is shown FIG. 1B. In this dual configuration, real projector 11 of FIG. 1A is transformed into a virtual camera 11″, and real camera 15 of FIG. 1A is transformed into a virtual projector 15″. It is to be understood that virtual camera 11″ and virtual projector 15″ represent the dual counterparts of real projector 11 and real camera 15, respectively, and are not real devices themselves. That is, virtual camera 11″ is a mathematical representation of how a hypothetical camera (i.e. virtual camera 11″) would behave to capture a hypothetically projected image 13″, which is similar to real image 13 projected by real projector 11 of FIG. 1A. Similarly, virtual projector 15″ is a mathematical representation of how a hypothetical projector (i.e. virtual projector 15″) would behave to project hypothetical image 13″ that substantially matches real image 13, as captured by real camera 15 (of FIG. 1A). Thus, the positions of the real projector 11 and real camera 15 of FIG. 1A are interchanged in FIG. 1B as virtual camera 11″ and virtual projector 15″. However, the pixel resolution of the real devices carries forward to their counterpart, dual devices. Therefore, virtual camera 11″ has a virtual receptor pixel array 17″ consisting of s rows and r columns to match projector pixel array 17 of real projector 11. Similarly, virtual projector 15″ has a virtual projection pixel array 19″ consisting of u rows and υ columns to match receptor pixel array 19 of real camera 15.

Suppose a dual transport matrix T″ is the transport matrix in this dual configuration such that a virtual captured image C″ (as captured by virtual camera 11″) is related to a virtual projected image P″ (as projected by virtual projector 15″) as C″=TP″, then T″_(eg) would be the dual light transport coefficient between virtual projector pixel g″ and virtual receptor pixel e″.

Helmholtz reciprocity specifies that the pixel-to-pixel transport is equal in both directions (i.e. from real projector 11 to real camera 15, and from virtual projector 15″ to virtual camera 11″). That is, T″_(eg)=T_(ge), which means T″=T^(T), (i.e. dual matrix T″ is equivalent to the result of the mathematical transpose operation on real matrix T). Thus, given matrix T, one can use T^(T) to synthesize the images that would be acquired in the dual configuration.

OBJECTS OF THE INVENTION

It is an object of the present invention to simplify the application of dual photography.

It is another object of the present invention to simplify the use of dual photography in the calibration of a real projector in projector-camera systems.

It is a further object of the present invention to make use of the benefits of dual photography in the calibration of a real projector while avoiding the need of some of the more time consuming and complicated requirements of dual photography.

It is a further object of the present invention to reduce some of the more time consuming and complicated requirements of feature identification and image system calibration.

SUMMARY OF THE INVENTION

The above objects are met in a system that simplifies the generation of dual transport matrix T, simplifies the implementation of light transport coefficients in dual photography, and combines the ease of dual photography manipulation with the reduced calculation requirement of homography.

Application of dual photography is simplified by reducing the number of captured images needed to generate a light transport matrix T of (p×q) projection pixel array from (p×q) images to (p+q) images. Manipulation of the light transport matrix is simplified by replacing the use of a matrix T with an index associating each projection pixel to only non-zero light transport values. By eliminating the use of zero-valued light transport coefficients, the memory and processing requirements for implementing dual photography are greatly reduced. This dual photography technique is applied to the calibration of projector-camera systems.

A second method for calibrating projector-camera systems uses a known projection test pattern and a calibrated camera to associate projected spatial markers on a real image to a captured image. Since the real location of the projected spatial markers are determined by the calibrated camera, and the relative location of the spatial markers on the test pattern is also known (since it was projected by the projector to be calibrated), the projector may be calibrated using homography techniques to assure that its known relation between projected spatial markers agrees with the true location of the spatial markers on the projection scene, as determined by the calibrated camera.

In a first embodiment of the present invention, a method of generating light transport coefficients relating a digital projector to a digital camera is shown. In the case where digital projector has an array of projection pixels and the digital camera has an array of sensor pixel, the method includes: simultaneously activating a first group of projection pixels within the projector to project a first test pattern on a projection scene, any projection pixels not in said first test pattern being maintained dark; capturing a first image of the first test pattern on the projection scene; simultaneously activating a second group of projection pixels within the projector to project a second test pattern on the projection scene, any remaining projection pixels not in the second test pattern being maintained dark, wherein the first and second groups of projection pixels have only one projection pixel in common defining a target projection pixel; capturing a second image of said second test pattern on said projection scene; comparing image pixels of the first image to corresponding image pixels of the second image and retaining the darker of two compared image pixels, the retained image pixels constituting a composite image; and identifying all none-dark image pixels in the composite image, the none-dark image pixels defining non-zero light transport coefficients associated with the target projection pixel.

Following this method, it may be further clarified that the first test pattern and the second test pattern form a first pattern pair associated with the target projection pixel; each projection pixel within the digital projector has an associated unique pattern pair comprised of a respective first test pattern and respective second test pattern; and the light transport coefficients of any selected projection pixel within the digital projector is obtained by comparing captured images of the selected projection pixel's respective first test pattern and respective second pattern comprising its associated unique pattern pair; wherein the comparing of the unique pattern pair includes: comparing image pixels of the respective first image to corresponding image pixels of said respective second image and retaining the darker of two compared image pixels, the retained image pixels constituting a respective composite image; identifying all none-dark image pixels in the respective composite image, the none-dark image pixels defining the light transport coefficients associated with any selected projection pixel.

Further in this method, the projector preferably has an array of rows and columns of projection pixels; the respective first test pattern is comprised of the simultaneous activating of all projection pixels within one column of the array; and the respective second test pattern is comprised of the simultaneous activating of all projection pixels within one row of the array.

Alternatively, the method of generating light transport coefficients may further include identifying the light transport coefficients for a selected number of the target projection pixels; generating an index associating each of the selected number of target projection pixels to their correspondingly associated non-zero light transport coefficients; and storing only the non-zero light transport coefficients.

In still an alternate embodiment, the light transport coefficients relating the digital projector to the digital camera are suitable for use in dual photography.

The present invention may be applied to the calibrating a projector-camera system having at least one digital projector and one digital camera, wherein the method includes: identifying a projection scene; calibrating the camera to the projection scene; identifying non-zero light transport coefficients associating each projection pixel within the digital projector to the digital camera; creating an index associating each projection pixel to its corresponding non-zero light transport coefficients, and disregarding each zero light transport coefficients of each projection pixel; using non-zero light transport coefficients, as identified by the index, to generate a dual image of the projection scene from the point of view of the projector using a dual photography technique; and applying a camera calibration technique to the projector by treating said projector as a virtual camera and using the dual image as a test image for calibrating said projector.

Further in this method, the step of identifying non-zero light transport coefficient associating each projection pixel within the digital projector to the digital camera, may include: identifying a first set of first test patterns, each first test pattern including of a plurality of projection pixels; identifying a second set of second test patterns, each second test pattern including a plurality of projection pixels, wherein each projection pixel within the projector is uniquely identified by an intersection of a test pattern pair consisting of a corresponding first test pattern from said first set and a corresponding second test pattern from said second set; activating each of the first test patterns within the first set onto the projection scene, and capturing an image of the resultant projected image; activating each of the second test patterns within the second set onto the projection scene, and capturing an image of the resultant projected image; and for each projection pixel within the projector, comparing image pixels of the captured image of its corresponding projected first test pattern to corresponding image pixels of the captured image of its corresponding projected second test pattern, retaining the darker of the two compared image pixels, the retained image pixels constituting a composite image, and identifying all none-dark image pixels in said composite image, the none-dark image pixels defining non-zero light transport coefficients.

In an alternate embodiment, a method for calibrating a projector-camera system having at least one digital projector and one digital camera, includes: identifying a projection scene; calibrating the camera to the projection scene; projecting a known test pattern unto the projection scene, the known test pattern including a group of spatial markers with a known spatial relation to each other; using homography techniques to identify the relationship between the camera and the projection scene to identifying the real locations of said spatial markers on the projection scene; using homography techniques to relate the real locations of the spatial markers on the projection scene as determined using said camera, to the known spatial relation of the spatial markers on the known test pattern, calibrate the projector so as make the known spatial relation of the known test pattern match the real locations of said spatial markers on said projection scene.

Continuing in this method, the known test pattern may include not less than four spatial markers.

Also, the known spatial relation between the spatial markers is preferably limited to a two-dimensional projection.

Further in the present invention, an image created by a first projector is recreated using a second projector by relating the two projectors to a common point of view, preferably as viewed from a camera. A first transport matrix T₁ is captured to relate the first projector to the camera. A second transport matrix T₂ is then capture to relate the second projector to the camera. To have a first image p₁, as projected by the first projector, reproduced by the second projector, the second projector projects a distorted image defined as (T₂ ⁻¹)(T₁)p₁.

If the first projector is front projector and the second projector is an immersive projector, then the immersive projector can be used to recreate a front projected image from the point of view of an audience without worrying about the audience blocking the projection path of a front projector.

In the calculation, (T₂ ⁻¹)(T₁)p₁, the inverse of T₂ is an estimation defined by first creating an intermediate matrix {hacek over (T)} of equal size as T₂. If a column in T₂ is denoted as Tr and a corresponding column in {hacek over (T)} is denoted as {hacek over (T)}r, then the construction and population of {hacek over (T)} is defined as {hacek over (T)}r=Tr/(∥Tr∥)², and the inverse of T₂ is estimated as the transpose of {hacek over (T)}.

The present approach may also be used to recreate real images, or lighting effects, first created in a computer model. In this case, a transport matrix relating a projector and a camera, which defines a preferred viewing position, is created. A model, such as a computer model of the projector, camera, and a projection surface is then created. Within the model, a model image or model lighting effect is created, and a virtual image as taken by the model camera is then generated. The virtual image is then generated using the real projector and real camera using the created transport matrix.

In embodiments, a method for identifying features in an image includes generating an image comprising a set of features and color patterns. Each feature in the set of features is associated with a color pattern that identifies the feature. The image is formed from at least one display image comprising a plurality of features or feature components used to form features. At least a portion of the image is captured and at least some of the features from the set of features in the at least a portion of the image are identified using color patterns.

In embodiments, a method for identifying features in an image includes using a projector to display a first display image comprising a first set of lines wherein each line has both a distinct color and a distinct position within the first set of lines and using a camera to capture a first image of at least a portion of the first display image. The projector is used to display a second display image comprising a second set of lines that are non-parallel to the first set of lines and wherein each line has both a distinct color and a distinct position within the second set of lines and the camera is used to capture a second image of at least a portion of the second display image. An image is generated from the intersections of the first set of lines and the second set of lines, wherein the intersections form features. Using color patterns in the image, at least some of the features can be identified.

In embodiments, a method for identifying features in an image includes generating an image comprising a set of features and color patterns. Each feature in the set of features is associated with a color pattern that identifies the feature. Because color information is used to help identify features, the image may be formed from at least one display image that has a plurality of features or feature components used to form features. At least a portion of the generated image is captured and at least some of the features from the set of features in the at least a portion of the image are identified using color patterns.

In embodiments, having identified at least some of the features in the captured image, an imaging system comprising a projector used to form the image and a camera used to capture the at least a portion of the image can be calibrated. In embodiments, one or more the methods previously described may be used for calibration, including using a homography transform for calibrating between the projector used to form the image and the camera used to capture the at least a portion of the image.

In embodiments, a method for calibrating an image system comprising a projector and a camera includes displaying, using the projector, at least one display image. The at least one displayed image is configured to form a calibration image comprising a set of features and color patterns. Each feature in the set of features is associated with a color pattern that identifies the feature. Also, the at least one displayed image is configured such that a plurality of features or feature components used to form features are displayed in a single image. Using the camera, at least a portion of the calibration image is captured and at least some of the features in the at least a portion of the calibration image are identified using color patterns. Having identified at least some of the features, those features can be used to calibrate the image system.

In embodiments, calibrating the system includes using position information of the identified features and a homography transform for geometric calibration between the projector and the camera used to capture the at least a portion of the image.

In embodiments, the at least one display image is a single image that forms the calibration image.

In embodiments, a method for calibrating an image system includes using a projector to display a first display image comprising a first set of lines wherein each line has a distinct position within the first set of lines, and a camera is used to capture a first image of at least a portion of the first display image. The projector is use to display a second display image comprising a second set of lines that are non-parallel to the first set of lines and wherein each line has a distinct position within the second set of lines, and the camera is used to a second image of at least a portion of the second display image. The projector is used to display a color image comprising color patterns, wherein a color pattern is uniquely associated with a feature position, and the camera captures a third image of at least a portion of the color image. It shall be noted that the order of the display and capturing is not critical to the present invention. A calibration image or at least a portion of a calibration image can be generated comprising intersections of the first set of lines and the second set of lines from the first and second images and the third image, wherein the intersections form features that are individually identifiable by the color patterns in the third image.

In embodiments, a computing system for calibrating an image system, which comprises a projector and a camera, includes an image generator that causes the projector to display at least one display image. The at least one displayed image is configured to form a calibration image comprising a set of features and color patterns. Each feature in the set of features is associated with a color pattern that identifies the feature. The at least one displayed image is also configured such that more than one features or feature components is displayed in at least one of the at least one display image. The system also includes a captured image processor that receives from the camera at least a portion of the calibration image and identifies at least some of the features in the at least a portion of the calibration image using color patterns. And, the system also includes a calibration processor communicatively coupled to the image generator and the captured image processor that uses position information of the identified features to calibrate the image system.

In embodiments, calibrating the system includes using position information of the identified features and a homography transform for geometric calibration between the projector and the camera of the image system.

In embodiments, one or more of the methods may be embodiment on a computer readable medium comprising one or more sequences of instructions which, when executed by one or more processors, causes the one or more processors to perform the one or more methods.

Other objects and attainments together with a fuller understanding of the invention will become apparent and appreciated by referring to the following description and claims taken in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

Reference will be made to embodiments of the invention, examples of which may be illustrated in the accompanying figures. These figures are intended to be illustrative, not limiting. Although the invention is generally described in the context of these embodiments, it should be understood that it is not intended to limit the scope of the invention to these particular embodiments. In the drawings wherein like reference symbols refer to like parts.

FIGS. 1A and 1B show a prior art setup for implementation of dual photography.

FIGS. 2A and 2B show a setup for dual photography in accord with the present invention.

FIG. 3A is an illustrative example of a projected footprint on a light sensor array within a digital camera resulting from activation of a single projection pixel in a projector.

FIG. 3B is an illustrative example of a column of light transfer coefficients within a matrix T reflecting the example of FIG. 3A.

FIGS. 4A and 4B show two examples of two columns of projection pixels simultaneously projected onto a scene having a checker board pattern.

FIGS. 5A and 5B show two examples of two rows of projection pixels simultaneously projected onto a scene having a checker board pattern.

FIG. 6A shows a generated light beam footprint resultant from a single projection pixel, as created by combining the images of FIGS. 4A and 5A.

FIG. 6B shows a generated light beam footprint resultant from a single projection pixel, as created by combining the images of FIGS. 4B and 5B.

FIG. 7 is first example of an index associating projection pixels to non-zero valued light transport coefficients, as determined by a light beam footprint as generated in FIG. 6A or 6B.

FIG. 8 is second example of an index associating projection pixels to non-zero valued light transport coefficients, as determined by a light beam footprint as generated in FIG. 6A or 6B.

FIG. 9A shows a real captured image taken by a real camera.

FIG. 9B shows a dual captured image, as seen by a real projector, as generated using the method of the present invention.

FIG. 10 is resultant example of application of the use of homography to calibrate a projector.

FIG. 11 is an exemplary projection setup for using dual photography to create an immersive display system.

FIG. 12 shows an immersive projector P2 being used to simulate a front projector P1″.

FIG. 13 is an example of an image generated using the virtual projector implementation of FIG. 12.

FIG. 14A shows the right side of an image projected by a real front projector.

FIG. 14B shows the corresponding left side of the image shown in FIG. 14A, but in FIG. 14B the left side of the image is projected by an immersive projector.

FIG. 14C shows the right side image of FIG. 14A joined to the left side image of FIG. 14B.

FIGS. 15A and 15B two additional examples of a left side image generated by an immersion projector joined to the right side image generated by a front projector.

FIGS. 16A to 16C show an alternate application of the present invention to recreate in a real room a virtual image created in a virtual model room.

FIG. 17 is an example of the technique of FIGS. 16A-16C to project an image bigger than the projection space in a real room without image distortion.

FIG. 18 is an exemplary projection system in accord with the present invention in its minimal form.

FIG. 19 shows a prototype based on the design of FIG. 18.

FIG. 20 is an alternate view of the setup of FIG. 19.

FIG. 21A shows, under ambient lighting, a room with the projection system of FIGS. 19 and 20 installed.

FIG. 21B shows the room of FIG. 21A under immersive projection lighting.

FIG. 22 is an alternative design for ceiling-mounted operation.

FIG. 23 is an alternate configurations of the present invention.

FIG. 24 is still another configuration of the present invention.

FIG. 25A shows a design that uses a single curved mirror 125 and multiple projector-camera pairs 145.

FIG. 25B shows a design that uses a single mirror pyramid 151 and multiple projector-camera pairs 145.

FIG. 26 shows how multiple large FOV projectors 153 a and 153 b can be used to achieve an even larger overall projection FOV.

FIG. 27 depicts an image system according to embodiments of the invention.

FIG. 28A depicts a system for calibrating an image system by identifying features using color information according to embodiments of the invention.

FIG. 28B depicts a computing system according to embodiments of the invention.

FIG. 29 illustrates a color-based method for calibrating an image system according to embodiments of the invention.

FIG. 30 depicts a set of colors that may be used in a color-based method for calibrating an image system according to embodiments of the invention.

FIG. 31 illustrates a color lines method for calibrating an image system according to embodiments of the invention.

FIGS. 32A and 32B illustrate images of color lines that may be used for calibrating an image system according to embodiments of the invention.

FIG. 33 illustrates a color grid method for calibrating an image system according to embodiments of the invention.

FIGS. 34A and 34B illustrate examples of color grids that may be used for calibrating an image system according to embodiments of the invention.

FIG. 35 illustrates an alternative color grid method for calibrating an image system according to embodiments of the invention.

FIGS. 36A and 36B illustrate examples of lines, which may be monochromatic lines that may be used in conjunction with a color grid for calibrating an image system according to embodiments of the invention.

FIG. 36C illustrates features obtained from the intersection of the lines displayed in FIGS. 36A and 36B that may be used in conjunction with a color grid for calibrating an image system according to embodiments of the invention.

FIG. 36D illustrates an overlay of the features obtained from the intersection of the lines displayed in FIGS. 36A and 36B and a color grid that may be used for calibrating an image system according to embodiments of the invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

If projectors could capture images, then there would be no need to differentiate between real projectors and real cameras, and projector-camera systems could be treated like multi-camera systems. Subsequently, the standard camera calibration techniques described above could be used to calibrate projector-camera systems. In other words, if the projector could be treated as a pseudo-camera, then it could be calibrated along with a real camera in a manner similar to the camera calibration stage of the multi-camera system described above, and the “bootstrapping” projector calibration stage might be eliminated.

Consequently, whereas in the “bootstrapping” projection calibration technique, a projector would project known feature points to be captured by a stereo pair of pre-calibrated cameras onto a known projection surface to determine the contours of the projection surface, a dual photography technique could eliminates the need for the stereo pair pre-calibrated cameras.

With reference to FIG. 2A, an imaging setup in accord with the present invention may include a real projector 21 and a real camera 25. Real projector 21 is preferably a digital projector and has an imaging element including an imaging projection array 27, consisting of p rows and q columns of individual imaging projection elements (i.e. projection pixels) j. Imaging projection array 27 is internal to real projector 21, and is shown in thin lines in a dotted box in FIG. 2A. Real projector 21 is preferably of the liquid crystal display (LCD) type, Digital Light Processing (DLP) type, liquid crystal on silicon (LCOS) type, or other digital projection technology type.

Preferably, real camera 25 is a digital camera having an image sensor including an image sensor array 29, consisting of m rows by n columns of individual image sensor elements (i.e. sensor pixels) i. For simplicity, image sensor array 29 is shown on real camera 25, but it is to be understood that image sensor array 29 is internal to real camera 25.

This physical setup using a real projector 21 and real camera 25 is preferably called the ‘primal’ setup. Light rays emitted from real projector 21 bounce off a projection surface (i.e. display environment or scene) 23, which may have an irregular or flat shape, and some of the light rays eventually reach the image sensor within real camera 25. In general, each ray of light is dispersed, reflected, and refracted in the scene and hits the camera's image sensor at a number of different locations throughout image sensor array 29. Thus, when a light ray emitted from an imaging projection pixel j in real projector 21 reaches real camera 25, the projected light ray forms an m-by-n image on the camera's image sensor array 29, with each individual sensor pixel i within image sensor array 29 receiving a certain amount of light. Consequently, each light ray emitted from each individual projection pixel j of real projector 21 generates (m×n) [i.e. m-multiplied-by-n] individual light transport coefficients, one for each sensor pixel i within real camera 25. Consequently, the entire light transport matrix T will consist of (p×q) [i.e. p multiplied-by q] columns (one for each individual projection pixel j within real projector 21) and (m×n) rows (one for each individual image sensor pixel i within real camera 25).

If a full, real image (i.e. an image projected using the entire imaging projection array 27 of real projector 21) is represented as a “p by q vector” [i.e. an (p×q) vector], Rprjct′, and the corresponding full, real image captured by the image sensor array 29 of real camera 25 is represented as an “m by n vector” [i.e. an (m×n) vector], Rcptr′, then the light transport between real projector 21 and real camera 25 can be written as Rcptr′=TRprjct′ where T is called the light transport matrix. It is to be understood that since each projection pixel j results in a light ray that is scattered across the entire image sensor array 29 of real camera 25, each sensor pixel i within the image sensor array 29 will have a differently valued light transport coefficient for each projection pixel j. Therefore, each projection pixel j will result in a column of (m×n) individual light transport coefficients, one for each sensor pixel i of the camera's image sensor array 29. Since there are (p×q) projection pixels j in real projector 21, light transport matrix T will have (p×q) columns [one for each projection pixel j] and (m×n) rows [one for each sensor pixel i] of individual light transport coefficients. Thus, light transport matrix T will necessarily be huge, consisting of (p×q×m×n) individual light transport coefficient values.

With reference to FIG. 2B, a “dual” setup is one where real projector 21 is replaced by a virtual camera 21″ having a virtual image sensor array 27″ of equal size as real imaging projection array 27. Thus, virtual camera 21″ has an image sensor array 27″ comprised of p-rows by q-columns of image sensor elements j″.

Similarly, in the “dual” setup, real camera 25 is replaced by a virtual projector 25″ having an imaging projection array 29″ of equal size as the image sensor 29. Therefore, virtual projector 25″ has an imaging projection array comprised of m-rows by n-columns of image projection elements i″.

In this case, a virtual image 23″ (as projected by virtual projector 25″) would be represented by a (m×n) vector, Vprjct″. Similarly, a virtual captured image captured by virtual camera 21″ could be represented by an (p×q) vector, Vcptr″. By the principle of Helmholtz reciprocity, the light transport is equal in both directions (i.e. from real projector 21 to real camera 25, and from virtual projector 25″ to virtual camera 21″). Therefore, one can write the dual light transport matrix, T″, for the dual setup (i.e. the duality transformation setup) as Vcptr″=T ^(T) Vprjct″ Remarkably, the duality transformation between the real light transport matrix, T, and its dual light transport matrix, T^(T), is the transpose matrix operation known in the art of matrix computation.

As it is known in the art of matrix computation, the transpose operation of a general [x×y] matrix A is denoted by A^(T) and is defined by an [y×x] matrix whose first column is the first row of matrix A, and whose second column is the second row of matrix A, whose third column is the third row of matrix A, and so on. As is readily evident, this matrix operation simply flips the original matrix A about its first element, such that its first element (i.e. at position (1,1)) remains unchanged and the bottom of the first column becomes the end of the first row. Consequently, if one can captures, or otherwise determines, the real light transport matrix T for a primal setup, then the dual light transport matrix T^(T) for the dual setup is readily computable by flipping the real light transport matrix T, as described.

As is explained above, real light transport matrix T holds the individual light transport coefficients corresponding between each individual projection pixel j of real projector 21 and all the individual sensor pixels i of real camera 25. Therefore, a determination of each individual light transport coefficient corresponding between an individual projection pixel j and all sensor pixels i should avoid light ray contributions from other projection pixels j in imaging projection array 27.

To accomplish this, one may first consider a full, real projection image [i.e. a (p×q) projection vector Rprjct′] that has zero values for all elements in imaging projection array 27 except for the j^(th) pixel (i.e. the projection pixel under-test for which the light transport coefficients are to be determined), which is set to value 1 (i.e. is turned ON). Then the captured, real image vector [i.e. the (m×n) vector Rcptr′], which is defined as Rcptr′=TRprjct′ is exactly the j^(th) column of matrix T. It is to be understood that since the resolution of real projector 21 is likely to be lower than the resolution of real camera 25, the single light beam emitted from the j^(th) projection pixel under test may fall across several sensor pixels i of image sensor array 29.

An example of this is shown in FIG. 3A where a partial view of image sensor array 29 shows individually sensors pixels i numbered horizontally from 1 to n on the first row, continuing with (n+1) to (2n) on the second row, and (2n+1) to (3n) on the third row, and so on. Following this sequence, it is to be understood that sensor pixels i along the bottom-most row would be numbered from (m−1)n+1 to (mn).

A light ray footprint 50 of a single light ray from a single projection pixel j impacting sensor array 29 is denoted as a circle. For illustration purposes, those sensor pixels i not within light ray footprint 50 [i.e. those sensor pixels i not hit by the single projection beam emitted from the j^(th) projection pixel] are shown as deeply darken, those pixels i partly covered by footprint 50 are shown as lightly darken, and those pixels i completely within footprint 50 are shown as having no darkening. As it is known in the art, each sensor pixel i that is at least partially covered by light ray footprint 50 will register a light intensity value proportional to the amount of light it receives. This light intensity value may be assigned as the light transfer coefficient for that individual sensor pixel i. Alternatively, the light transport coefficient of each sensor pixel i may be made proportional to the light intensity value registered by the individual sensor pixel i. Nonetheless, those sensor pixels i that are not directly hit by the projection beam from the j^(th) projection pixel will have a value of zero, and thus their corresponding light transport coefficient will likewise have a value of zero.

With reference to FIG. 3B, an example of an image capture vector Rcptr′_j [or j^(th) column of matrix T], as might correspond to the footprint example of FIG. 3A is shown. This j^(th) column of matrix T is illustratively shown as an numbered sequence of vertical light transport coefficients, each corresponding to the numbered sensor pixels i of FIG. 3A. The numerical sequence of capture vector Rcptr′_j preferably follows the horizontally numbered sequence of individual sensor pixels i in image sensor array 29 shown in FIG. 3A. As shown, only those elements in captured imaged vector Rcptr′_j that correspond to sensor pixels covered by light ray footprint 50 have non-zero, i.e “NZ”, values for light transport coefficients. It is to be understood that “NZ” represents any non-zero light coefficient value, and that this value would be related to the amount of light intensity received by the corresponding sensor pixel i. Since light ray footprint 50 spans several rows of image sensor array 29, and each row is sequentially listed in captured image vector Rcptr′_j, several long series of zero valued light transport coefficients may be interspersed between few non-zero, NZ, valued light transport coefficients.

Nonetheless, this method of acquiring a column of light transport coefficients for matrix T for a given j^(th) projection pixel, suggests that a systematic method for capturing the entire matrix T is to sequentially turn ON each projection pixel j of real projector 21 (one projection pixel at a time), and to capture its corresponding real image, Rcptr′_j, with real camera 25. When all p×q projection pixels j have been sequentially turned ON, and their corresponding real images Rcptr′_j have been captured, all the captured image vectors Rcptr′_(1-to-(p×q)) are assembled into matrix T. Each captured image vector Rcptr′_j constitutes a column of light transport coefficient entries in matrix T. This results in a matrix T having (p×q) columns and (m×n) rows of individual light transport coefficients.

This straight forward, and systematic process for determining matrix T, however, is obviously a time-consuming process requiring (p×q) image projection-and-capture steps. Furthermore, the resultant light transport matrix T is very large, consisting of (p×q×m×n) elements. Because of the size of matrix T, computing a dual image is an extremely computation intensive operation requiring matrix multiplication between dual light transport matrix, T^(T), (which has the same number of elements as matrix T), and virtually projected image vector Vprjct″ (which is a long vector having (m×n) elements).

In the past, a scheme for determining matrix T by adaptively turning on multiple pixels of a real projector at a time (such as projector 11 of FIG. 1A), has been suggested to speed up the process of determining matrix T. However, this scheme is complex, and requires that a projection image be divided into multiple zones, and that one pixel from each zone be selected for simultaneous lighting while assuring that simultaneously lit projection pixels are maintained distant enough from each other to eliminate light interference between them. Additionally, this scheme is scene-dependent (i.e. dependent upon the display environment) since it is necessary to assure that interference between individually projected light rays is minimized. Thus, creation of the multiple zones is dependent on the display environment. Therefore, this scheme requires much physical inspection of a projection area and manual setup. Consequently, this scheme does not easily lend itself to general use.

A feature of the present invention proposes a method of reducing the number of procedural steps in the determination of real light transport matrix T. That is, instead of requiring (p×q) image projection-and-capture steps (and the storing of p×q captured images), the presently proposed method captures only “p plus q” [i.e. (p+q)] images.

The present method is based on the following assumptions: for most projector-camera display applications, any two distinct light rays b and c emitted from real projector 21 will typically hit image sensor array 29 in real camera 25 at distinct parts. That is, the overlap in the sensor pixels i hit by light from each of the projected light rays is assumed to be small, and may be ignored.

Of course, this is not true in general. For example, if the scene (i.e. display environment or projection surface/area) consists of a glass of milk (or other light diffusing material), the projected light rays will be diffused by the milk and there will be significant light overlap between the different light rays at the camera's image sensor array. However in a display setup designed to ensure high resolution projections it is virtually guaranteed that each projected pixel will be distinct from the next. That is, in venues, or settings, where high resolution projections are desired, it is likely that the venue will be clear of light diffusing articles along the light path of a projected image.

Continuing with the present novel method for determining matrix T, it is first noted that individual light transport coefficient contributions from projector pixel j map to column j of the light transport matrix T. Assuming minimal overlap between projection pixels, it follows that a first set of projection pixels S1 within imaging projection array 27 [i.e. S1⊂{1, . . . , (p×q)}] maps to a corresponding set of columns (one per projection pixel) in light transport matrix T. Furthermore, it is assumed that the first set of projection pixels S1 includes target projection pixel j, i.e. the target pixel under test.

Let Rcptr_S1 be a first image captured by real camera 25 of a projected image created by the simultaneous activation of first set of projection pixels S1 .

Consider now a second set of projection pixels S2 who share only projection pixel j in common with the first set of projection pixels S1 , [i.e. S1∩S2={j}] Let Rcptr_S2 be a second image captured by real camera 25 of a projected image created by the simultaneous activation of the second set of projection pixels S2 . The light transport coefficients of the j^(th) column of light transport matrix T (which corresponds to the target projection pixel under test, i.e. correspondents to j) may be directly obtain from real captured images Rcptr_S1 and Rcptr_S2 by identifying the one light ray footprint (i.e. similar to light ray footprint 50 in FIG. 3A) they share in common. This common light ray footprint would correspond to a light ray emitted from target projection pixel j, which is the only lit projection pixel that first set S1 and second set S2 share in common.

A method of identifying this common light ray footprint is to conduct a pixel-by-pixel comparison of both captured images Rcptr_S1 and Rcptr_S2 . For example, in first captured image Rcptr_S1 only sensor pixels within light ray footprints, each corresponding to the simultaneous lighting of the first set of projection pixels S1 , will have non-zero (NZ) light intensity values, and all other pixels in captured image Rcptr_S1 will have zero values, i.e. will be comparatively dark. Similarly in second captured image Rcptr_S2 , only sensor pixels within light ray footprints corresponding to second set of projection pixels S2 have non-zero (NZ) light intensity values, and all other pixels will have zero (or dark) values (i.e. below a predefined threshold value). Since the two sets S1 and S2 share only the target projection pixel, j, in common, a direct comparison of both captured images will quickly identify the sensor values corresponding to projection pixel j by identifying the only non-zero region (i.e. non-dark region) common to both Rcptr_S1 and Rcptr_S2 . Stated differently, the intersection of the lit regions (i.e. light ray footprints) of Rcptr_S1 and Rcptr_S2 is identified, and this intersection corresponds to the target projection pixel, j.

A method of accomplishing this is to conduct a pixel-by-pixel comparison of both captured images Rcptr_S1 and Rcptr_S2 , and retaining only the darker of the two compared pixels. This process may be expressed as: Tj≈MIN(Rcptr _(—) S1,Rcptr _(—) S2) where Tj is the j^(th) column of matrix T, and “MIN” indicates that the lower valued sensor pixel (i.e. the darker pixel having a lower captured light intensity value) in Rcptr_S1 and Rcptr_S2 is retained, and the higher valued (i.e. brighter) sensor pixel is discarded. In this way, the only high intensity values that are retained correspond to a lit footprint common to both S1 and S2 .

Stated differently, since the contribution of each individual projection pixel j is mapped to distinct parts of the image sensor array 29 of real camera 25, there is a set of sensor pixels L among sensor pixels 1-to-(m×n), [i.e. L⊂{1, . . . , (m×n)}] common to captured image Rcptr_S1 and Rcptr_S2 that corresponds to the target projector pixel, j. It should again be noted that the target projection pixel, j, is the intersection of projection pixel sets S1 and S2 , (i.e. j is the only projection pixel common to both sets S1 and S2 ), such that S1∩S2={j} Therefore, among the captured image pixels (in both Rcptr_S1 and Rcptur_S2 ) that do not correspond to the target projection pixel, j, [i.e. those sensor pixels not in set L, i.e. ∉L], at least one of the compared, captured sensor pixels in either Rcptr_S1 or Rcptr_S2 will not have received light. Since sensor pixels receiving light will be brighter than sensor pixels not receiving light, the operation MIN(Rcptr_S1 , Rcptr_S2 ), provides an image where only pixels in set L [i.e. ∈L] are lit, which is a good approximation of Tj, i.e. the j^(th) column in matrix T.

This implies that if sets of adjacent projection pixels in imaging projection array 27 are lit in columns and in rows, and a first collection of captured images Rcptr_Sy are made for the lit columns of projection pixels and a second collection of captured images Rcptr_Sx are made for the lit rows of projection pixels, then the light coefficients for any individual projection pixel j may be obtained by comparing both collections and identifying the region L where a captured image of a lit column intersect a captured image of a lit row, the intersection corresponding to a light ray projected by activation of projection pixel j, alone.

Thus, a method of determining transport matrix T is to collect a set of images Rcptr_Sy_1, . . . , Rcptr_Sy_q, corresponding to q captured images of q lit columns of projection pixels, and construct a set Rcptr_Sx_1, . . . , Rcptr_Sx_p corresponding to p captured images of p lit rows of projection pixels. Then for all projection pixels j in the imaging projection array of real projector 21 [which consists of projection pixels 1 . . . (p×q)], there exists a pair of row and column image captured sets, Rcptr_Sy_a, Rcptr_Sx_b, such that the intersection region of sets Rcptr_Sy_a and Rcptr_Sx_b correspond to a light ray footprint created by activation of target projection pixel, j. Therefore, one needs to construct sets of projection images Rprjct _(—) Sy _(—)1, . . . , Rprjct _(—) Sy _(—) q and Rprjct _(—) Sx _(—)1, . . . , Rprjct _(—) Sx _(—) p where each image Rprjct_Sy_1, . . . , Rprjct_Sy_q is paired with any of images Rprjct_Sx_1, . . . , Rprjct_Sx_p such that each set of paired images shares only one projection pixel in common, such that ∀j∈{1, . . . , (p×q)}∃Rprjct _(—) Sy _(—) a, Rprjct _(—) Sx _(—) b|Rprjct _(—) Sy _(—) a∩Rprjct _(—) Sx _(—) b={j} The above formula being interpreted to mean that for all projection pixels j in {1 . . . (p×q)} there exist a pair of projection images, each having a differently constructed pattern such that the intersection of the constructed patterns intersect at a single point corresponding to a common projection pixel. A basic example of such pairs of constructed patterns would be projected pairs of vertical light beams and horizontal light beams. In this case, the intersection of a captured image of a vertical line and a captured image of a horizontal line would include all the sensor pixels i that correspond to a target projection pixel, j, (i.e. pixels i that lie within a light ray footprint created by a light ray emitted from projection pixel j).

Therefore, any column Tj [where j=1, . . . , (p×q)] within transport matrix T can be synthesize from images Rcptr_Sy_1, . . . , Rcptr_Sy_q and Rcptr_Sx_(—)1, . . . , Rcptr_Sx_p.

A scheme that satisfies this property is to use pixel coordinates: let Rprjct_Sx_j be a first projected image such that only pixels with an x-coordinate equal to j are turned ON, and Rprjct_Sy_k be a second projected image such that only pixels with an y-coordinate equal to k are turned ON. Then MIN(Rprjct_Sx_j, Rprjct_Sy_k) gives an image of only coordinates (j, k) turned ON. This process can best be understood with reference to FIGS. 4A, 4B, 5A, 5B, 6A, and 6B.

In FIG. 4A, a scene, or display environment, consists of a flat surface 41 with a checker-board pattern. The relative location of each dark block 43 and each light block 45 that make up the checker-board pattern on flat surface 41 are known, and their known relative displacement may be used to calibrate real projector 21, as is more fully explained below.

Firstly, a bright vertical beam (or column of light rays), 47 _(—) k is projected onto surface 41 by real projector 21. In the present case, vertical beam 47 _(—) k is generated by turning ON all projection pixels within the imaging projection array 27 that have a y-coordinate equal to k. Real camera 25 then captures this image, Rcptr_Sy_k, as one example of a lit column of projection pixels.

In FIG. 4B, where all elements similar to those of FIG. 4A have similar reference characters, real projector 21 projects a second vertical beam 47 _(—) t onto surface 41. In this case, vertical beam 47 _(—) t is generated by turning ON all projection pixels having a y-coordinate equal to t. Real camera 25 then captures this image, Rcptr_Sy_t, as another example of a lit column of projection pixels. It is to be understood that real projector 21 would project a separate vertical beam for each of the q columns in it (p×q) imaging projection array 27, and real camera 25 would capture a separate image of each projected vertical beam.

With reference to FIG. 5A, all element similar to those of FIGS. 4A and 4B have similar reference characters and are described above. In the present case, real projector 21 projects a horizontal beam 49 _(—) j onto projection surface 41. Horizontal beam 49 _(—) j is generated by turning ON all projection pixels having an x-coordinate equal to j. Real camera 25 then captures this image, Rcptr_Sx_j, as one example of a lit row of projection pixels.

In FIG. 5B, real projector 21 projects a second horizontal beam 49 _(—) r onto surface 41. As before, horizontal beam 49 _(—) r is generated by turning ON all projection pixels having a x-coordinate equal to r. Real camera 25 then captures this image, Rcptr_Sx_r, as another example of a lit row of projection pixels. It is to be understood that real projector 21 would project a separate horizontal beam for each of the p rows in it (p×q) imaging projection array, and real camera 25 would capture a separate image of each projected horizontal beam.

If one now conducts a pixel-by-pixel comparison of captured image Rcptr_Sy_k from FIG. 4A and captured image Rcptr_Sx_j from FIG. 5A (or alternatively compares only their respective vertical and horizontal beams), using operation MIN(Rcptr_Sx_j, Rcptr_Sy_k) to retain only the darker of two compared image pixels and discarding the brighter of the two, one would generate an image 41′, as shown in FIG. 6A. All elements in FIGS. 6A and 6B similar to those of FIGS. 4A, 4B, 5A, and 5B have similar reference characters and described above, with the addition of an apostrophe.

Since most of the Rcptr_Sx_j image is the same as image Rcptr_Sy_k (i.e. they mostly consist of the plain checker board pattern on flat surface 41 with projected light rays), retaining the darker of two compared pixels does not change the majority of the resultant image. That is, if two compared pixels are relatively the same, then electing either pixel over the other does not much affect the resultant image 41′. However, when a pixel on vertical light beam 47 _(—) k in image Rcptr_Sy_k is compared with a corresponding pixel in image Rcptr_Sx_j that does not lie on horizontal light beam 49 _(—) j, then retaining the darker of the two image pixels will discard the bright, beam pixel from image Rcptr_Sy_k, and retain the plain pixel from image Rcptr_Sx_j, which shows no light beam. Therefore, vertical light beam 47 _(—) k is eliminated from generated image 41′. Similarly, when a pixel on horizontal light beam 49 _(—) j in image Rcptr_Sx_j is compared with a corresponding pixel in image Rcptr_Sy_k that does not lie on vertical light beam 47 _(—) k, then retaining the darker of the two pixels will discard the bright, beam pixel from image Rcptr_Sx_j, and retain the plane pixel from image Rcptr_Sy_k, which shows no light beam. Consequently, horizontal light beam 49 _(—) j is also eliminated from generated image 41′. However, within the region where vertical light beam 47 _(—) k intersects horizontal light beam 49 _(—) j, both compared image pixels are bright beam pixels showing an impact by a light ray. Comparison of these two image pixels within this intersection region will result in either of the two bright beam pixels being selected for image 41′. As a result, image 41′ will show a brightly lit region 53 corresponding to a projected light ray emitted from coordinates (j,k) of imaging projection array 27. Thus, the light transport coefficients for the projection pixel having coordinates (j,k) can be extracted from generated image 53 without having to have physically captured an image of a light ray projected from the projection pixel at (j,k).

A second example is shown in FIG. 6B, where the combination of real captured images corresponding to FIGS. 4B and 5B (which would respectively correspond to real captured images Rcptr_Sy_t and Rcptr_Sx_r following the above-described naming convention), results in a second brightly lit region 55 corresponding to a projected light ray emitted from coordinates (r,t) of image projection array 27.

A similar process may be followed to identify the light transport coefficients of every projection pixel in imaging projection array 27 without having to individually turn ON and project each projection pixel j, one-at-a-time. This method of generating an image of a hypothetically, singularly activated projection pixel to obtain the projection pixel's light transport coefficients requires only (p+q) captured images, one for each row and column of projection pixels in imaging projection array 27 of real projector 21. Furthermore, once all the pixel projection locations have been identified, the (p+q) captured images may be discarded, and all that needs to be saved is an index and corresponding footprint information.

An example of this approach is shown in FIG. 7, where an index of real imaging projection array 27 is represented as a partial array with circles 1, 2, 3, . . . (q+1) . . . (2q+1) . . . etc. representing individual projection pixels j. Similarly, real image sensor array 29 is represented as a partial array with circles F1, F2, F3, . . . F(q+1), . . . etc. representing the footprint information corresponding to individually activated projection pixels j. In the present example, footprints F1, F2, F3, . . . F(q+1), . . . etc. respectively correspond to projection pixels 1, 2, 3, . . . (q+1), etc. For illustrative purposes, the partial array representing sensor array 29 is made more denser than the array representing imaging projection array 27 in order to better illustrate that pixel density (i.e. resolution) of real camera 25 may be much greater than the resolution of real projector 21, and thus a light beam emitted from a single projection pixel j may create a light footprint spanning several sensor pixels i.

Following this approach, only two sets of information need to be stored. A first set of information corresponds to an index of projection pixels and a second set of information corresponds to sensor pixel footprint information corresponding to each projection pixel. In other words, zero coefficients need not be stored, which greatly reduces the memory requirements.

A second example of organizing this information is shown in FIG. 8, where an index 61 of projection pixel is shown to point, or correspond, to a groups 63 of grayscale (i.e. non-zero) sensor pixel information (i.e. corresponding to a resultant light beam footprint).

Having shown how to reduce the number of images that need to be captured and stored to generate the needed light transport coefficients in preparation for dual photography, the following now addresses some of the practical difficulties of implementing dual photography. As is explained above, a light transport matrix T can be very large, and its use (or the use of its transpose, the dual light transport matrix T^(T)) requires large amounts of active memory (for example, DRAM) and excessive computational processing power/time. Therefore, general use of the dual image has heretofore not been practical.

To efficiently construct a dual image, one first notes that Vcptr″=T ^(T) Vprjct″ Since the virtual image sensor array 27″ in the virtual camera 21″ corresponds in actuality to real imaging projection array 27 of real projector 21, it is convenient to use the same index j to denote any virtual sensor pixel within Vcptr″ obtained by virtual camera 21″. Therefore, a relationship between each virtual projection pixel j in a virtual captured image versus a corresponding row of elements in T^(T) may be denoted as Vcptr″(j)=T ^(T) _(j) Vprjct″ where T^(T) _(j) refers to the j^(th) row in T^(T).

As is explained above, T^(T) is the transpose of light transport matrix T (i.e. matrix T turned on its diagonal), and the values of row T^(T) _(j) (where j is any value from 1 to (p×q)) therefore correspond to the j^(th) column of matrix T (i.e. T_(COL) _(—) j). Since each column of T has (m×n) elements (i.e. equivalent to the pixel resolution of real camera 25), this would appears to be a very large number of elements. However, recalling that in the present implementation, only a limited number of elements in each column of matrix T are non-zero (i.e. only those corresponding to camera sensor pixels i upon which shone the intersection of a vertical and horizontal light beam, i.e. a light footprint), it is self apparent that only a few of the (m×n) elements within in each column j of matrix T (and subsequently in each row T^(T) _(j)) are non-zero. Therefore, it is not necessary to examine all elements in T_(COL) _(—) j when computing Vcptr″ (j). Indeed, as is explained above, it is preferred that only a single index 61 showing all individually projected pixels, and their corresponding light beam footprint information 63 be stored, as illustrated in FIGS. 8 and 7.

As is shown in FIG. 8, index 61 associates a listing of grayscale entries 63 for each projected light beam (from an individual projection pixel j of real projector 21). It is to be understood that each group of grayscale entries 63 corresponds to the non-zero entries within each column of T, and only these non-zero grayscale values need to be examined during each matrix operation of a column of T. Thus, the number of calculations needed to for determining each value of Vcptr″ (j)=T^(T) _(j)Vprjct″ is greatly reduced.

In other words, this subset of elements, SVcptr″ (G), within each column G of matrix T that needs to be processed may be defined as T_(COL) _(—) G(a), where a is an index for any virtually captured pixel [i.e. a∈{1, . . . , (p×q)}] in Vcptr″(G). Therefore, for each SVcptr″(G), one can define the set of elements to be examined as: SVcptr″(G)={a|∀

∈{1, . . . , (p×q)}T _(COL) _(—) G(a)≧T T _(COL) _(—)

(a)} Since in general ∥SVcptr″(G)∥<<(p×q), it takes significantly less time to compute:

${{Vcptr}^{\prime\prime}(G)} = {\sum\limits_{\alpha \in {{SVeptr}^{\prime\prime}{(G)}}}{{T_{G}^{T}(\alpha)}{{Vprjct}^{\prime\prime}(\alpha)}}}$ than to compute: Vcptr″(j)=T ^(T) _(j) Vprjct″

An example of a dual image generated using this method is shown in FIGS. 9A and 9B. FIG. 9A shows a primal image, as projected by a real projector. FIG. 9B shows the resultant dual image computed by an implementation of present method. The dual image of FIG. 9B represents the image virtually captured by the real projector 27 (i.e. virtual camera 21″), or stated differently, the image as “seen” by the real projector 21.

The above discussion shows how to compute dual images efficiently from a reduced set of images, which saves image capture time as well as computation time. As is explained above, the real captured images and dual captured images can be used to calibrate both real camera 25 and real projector 21, respectively.

That is, since the images from the projector's view cannot be directly captured, a straightforward solution is to construct projector-view images (i.e. dual images) from corresponding camera images using dual photography techniques, and then to calibrate the projector using the constructed images. For example, after having taken several real images at different angles of a known object with carefully measured features, the real camera can be calibrated by using the known dimensions of the object to compensating for distortions in the captured images arising from the different angle views. The virtual images, as seen by the real projector, can then be generated from the same captured images using dual photography techniques, as described above, and the real projector may be calibrated in a manner analogous to the real camera.

A possible set-back associated with this straight forward method, however, is the difficulty in generating and manipulating the light transport matrix T, and operating on the large image vectors resulting from the large number of camera and projector image pixels. Although this labor-intensive and expensive process is mitigated substantially by using the dual photography method described above, for purposes of calibrating real projector 21 in a projector-camera system, such as that shown in FIG. 2A, Applicants have developed a novel method that avoids the need for generating a T matrix and create a full dual image, while still taking advantages of some of the benefits of using a dual image (i.e. an image as “seen” by real projector 21) to ease calibration of real projector 21.

The generation of the T matrix can be avoided altogether by noting that to calibrate the projector, one does not need to construct an entire dual image, but only needs to determine the location of the corners of the squares within the checker board pattern on flat surface 41 (see FIG. 4A), as seen by the projector. Applicants have adapted a homography-based method to achieve this goal, and thus avoid the generation and manipulation of the T matrix of fully-applied, dual photography. This alternate embodiment of the present invention directly computes the coordinates of the checker corner features across the projector-view images without requiring the construction the dual images and the detection of the corners from the constructed dual images.

In this novel method, the real camera 25 is first calibrated using any known method, such as the previously described method of capturing multiple images of a known pattern at different angles. More specifically, one may use the real camera to capture a number of images of a planar checker board at different poses and detect the checker corners across the captured images. It is to be understood that that a checker board is being used purely for illustrative purposes, and any scene may be captured as long as physical distances between key points are known. In the present case, since the physical distances between the checkered corners on the board have been previously measured, their captured images at different poses are sufficient to compute the internal and external camera parameters, as was described above.

Secondly, it is observed that projector images follow the so-called perspective projection model, which relates two (or more) views of a single scene as seen by two (or more) separated sources. That is, different viewing sources will “see” a different view (or image) of the same scene since the different sources are located at different angles to the scene. However, since there is only one real scene (irrespective of the number of views of the scene), one can generate a mathematical relationship between the different views that will associate any point on any one view to a corresponding real point on the scene (and thereby to all other views).

If one of these separated views is assumed to be a virtual image as “seen” by a real projector, while a second separated view is deemed to be a real image captured by a real camera, then the perspective projection model (which relates the two views to the common, real scene) would permit one to extract from the captured real image some information relating to the virtual image, without generating a full dual image.

Using this approach, Applicants have devised a method of extracting sufficient information for calibrating a real projector without requiring a full dual image. Thus, although no dual image has been created, one can still apply a camera calibration technique to a projector, albeit in a round-about way.

Under the perspective projection model, the relationship between two image projections of a planar object from different views is a simple linear projective transformation or homography. This transformation relates the coordinates of any point on the planar object (i.e. a homogeneous coordinate) to the coordinates of a corresponding point on a specific view of the planar object. In the present embodiment, the projector-view image of the planar checker board is a homography of the corresponding camera image. Specifically, for any point P on the checker board, its homogeneous coordinate in the projector-view image Up=(up, υp, 1) and the coordinate in the camera image Uc=(uc, υc, 1) satisfy the following equation, Up=λHUc where λ is a scalar and H is the 3×3 homography transformation matrix (as known in the art) of which the bottom right entry is set to 1. The pair of corresponding coordinates provide 3 linear equations, where one of them determines the scalar and the other two are used to determine H, the homography transformation matrix. Since there are 8 unknown entries in 3×3 matrix H, given the correspondence of between N coordinate points (where N≧4) on the checker board, the homography between the projector-view image and the camera image can be recovered by solving the 2N linear equations. The greater the number of N, the lower the error relating coordinate points between the projector-view and the camera image.

To obtain the corresponding coordinates, 10 white points are preferably projected on the checker board pattern, an image of the checker board with the projected white points is captured using a real camera, such as real camera 25, and the coordinates of the 10 points in the camera image are computed. In the present process, it is only necessary that the ten points be distinguished during the computation of their corresponding coordinates in the captured image. This may be achieved by projecting the ten points sequentially, and determining their corresponding coordinates, in turn. Alternatively, differently colored points may be projected simultaneously, and the different points identified by color.

Since the projector preferably projected the points in a known relation to each other, the coordinates of the points in the projected image are known. This results in 10 pairs of corresponding coordinates, one set as captured by the real camera and a second set as projected by the real projector. Once the homography is recovered, the coordinates of the checker board corners detected in the camera images can be directly transformed to compute the corresponding corner coordinates in the projector-view images. The projector parameters can then be calibrated using a camera calibration method, such as the one described above.

An example of this approach is shown in FIG. 10. In FIG. 10, the feature capture results are as follows. The circles, or dots, not shown in outline (for example dots 81) were used for estimating homography, while the outlined circles, or dots, (for example dots 83) are the corner point features. As can be seen, the outlined dots 83 are on the actual corners, indicating that the projector coordinates for each detected corner has been correctly captured.

A still alternate embodiment may be obtained by taking another look at the setup of FIG. 2A. As is explained above, a real captured image (Rcptr′) is related to a real projected image Rprjct′ by the light transport matrix T. As is also explained above in reference to FIG. 2B, a virtual captured image Vcptr″, as captured by projector 21′ (illustratively shown as a virtual camera 21″) if a projector could capture images, is related to a virtual projected image Vprjct″, as projected by camera 25′ (illustratively shown as a virtual projector 25″) if a camera could project images, by the transpose of light transport matrix T (i.e. T^(T)). That is, Vcptr″=T^(T) Vprjct″. However, there are situations where using the inverse of transport matrix T, instead of its transpose, may be useful. For example, in article “A Theory of Inverse Light Transport” by Seitz et al. (IEEE International Conference on Communications, ICC V05), hereby incorporated in its entirety by reference, the inverse of the light transport matrix T was used in to analyze the way light bounces in arbitrary scenes. In this approach, a scene is decomposed into a sum of η-bounce images, where each image records the contribution of light that bounces η times before reaching a camera. Using a matrix of “impulse images”, each η-bounce image is computed to infer how light propagates through the scene.

Another example of how the inverse of transport matrix T may be used can be understood by reconsidering the projector-camera system of FIG. 2A, where Rcptr′=T Rprjct′. If one has an image (such as Rcptr′, for example) but wants to infer an unknown projection source that can produce the image, one can compute Rprjct′=T⁻¹ Rcptr′. Note that T^(T) Rprjct″ shows virtually captured image Vcptr″, which is a different view of Rcptr′, but does not show the real source of projected image Rprjct′.

T⁻¹ is however harder to compute than the transpose, requiring much more computational resources. Worse, it is not always possible to find the inverse of an arbitrary matrix. That is, some matrixes may not have an inverse. Fortunately, Applicants have found that in many display settings an approximation to the inverse of transport matrix T can be computed. As it is known in the art, the identity matrix, or unit matrix, is denoted by I, and has the property that for a given matrix A, the following relationship holds: AI=IA=A If matrix A were a matrix of order m by n, then the pre-multiplicative identify matrix I would be of order m by m, while the post-multiplicative identify matrix I would be of order n by n.

The multiplicative inverse of a matrix is typically defined in terms of identify matrix I. A left multiplicative inverse of a matrix A is a matrix B such that BA=I, and a right multiplicative inverse of a matrix A is a matrix C such that AC=I. If the left and right multiplicative inverse of a matrix A are equal, then the left and right multiplicative inverse is simply called the “multiplicative inverse”, or inverse, and is denoted by A⁻¹.

In general, T⁻¹≠T^(T). However, as discussed above, for most projector-camera display applications, any two distinct light rays j and k emitted from a projector will typically hit a camera sensor at distinct parts, i.e., there is usually little overlap in the camera sensor pixels hit by light from each of the light rays j and k. Of course, this is not true in general, such as the example given above where the scene consists of a glass of milk, and the light rays are diffused by the milk resulting in significant overlap. However in a display setup designed to ensure high resolution projections it is virtually guaranteed that each pixel projected will be distinct from the next.

As is also explained above, each column of the transport matrix T is the projection image of one pixel from the projector. Thus all of the column entries have zero values except those corresponding to the camera pixels hit by the projector pixel. Under this display setting, and with the assumption that camera pixels corresponding to different projector pixels do not overlap with each other, the columns of T are orthogonal to each other, meaning that that the transpose of a specific column multiplied by any column other than itself will produce a zero result.

To determine the inverse of T, i.e. T⁻¹, it is beneficial to first note that AA⁻¹=I, and the identify matrix I is comprised of a matrix with numeral 1's along a diagonal from the top left corner (starting at point 0,0) to the bottom right corner (ending at point r,g) and zeros everywhere else. Therefore, to compute T⁻¹, a matrix {hacek over (T)} is constructed such that each column in {hacek over (T)} is comprised of normalized values of corresponding values in a corresponding column in T. That is, {hacek over (T)}r=Tr/(∥Tr∥)² , r=1, 2, 3, . . . , pq where {hacek over (T)}r is the r^(th) column of {hacek over (T)}. Since the operation ∥Tr∥ defines the square root of the sum of the squares of all values in column r of matrix T, the square of ∥Tr∥ is simply the sum of the squares of all the values in column r. That is,

${{Tr}} = \left. \sqrt{}\left. \overset{\omega}{\underset{ɛ}{\{}}{\sum\left( {Tr}_{ɛ} \right)^{2}} \right\} \right.$ By dividing each value entry in column r by the sum of the squares of all the values entries in column r, operation {Tr/(∥Tr∥)²} has the effect of normalizing the value entries in column r of matrix T. If one now takes the transpose of {hacek over (T)}r, i.e. flips it on its side such that the first column becomes the top row and the last column becomes the bottom row, the result will be rows of elements that are the normalized values of corresponding columns of elements in T. Therefore, for every column in T, one has the following result: ({hacek over (T)}r ^(T))×(Tr)=1 and ({hacek over (T)}r ^(T))×(Tg)=0, for r≠g In other words, multiplying a column of T with a corresponding row in {hacek over (T)}r^(T) always results in numeral 1, and as one multiplies all the columns in T with the corresponding row in {hacek over (T)}r^(T), one produces a matrix with numeral 1's along its diagonal, and one may place zeroes everywhere else.

Therefore, in the case of matrix T, where columns are orthogonal to each other, and given the specific construction of matrix {hacek over (T)}, one has that the transpose of {hacek over (T)} is equivalent to the inverse of T (i.e. {hacek over (T)}^(T)=T⁻¹), by definition, or at least {hacek over (T)} is a left multiplicative inverse of T. Therefore, Rprjct′={hacek over (T)}^(T) Rcptr′.

Note that only part of the projector pixels that actually hit the camera sensor can be recovered. For the projector pixels not hitting any of the camera pixels, the corresponding columns in T contain purely zeros and the above equation of {hacek over (T)}r=Tr/(∥Tr∥)² is undefined. In such cases it is preferred that the corresponding columns in {hacek over (T)} be set as zero columns. Thus {hacek over (T)}^(T) is the inverse of the part of T that covers the overlapping area of the field-of-views of the projector and the camera. It only recovers the projector pixels in Rprjct′ that fall in the overlapping area and blacks out the other pixels.

Having defined a method for approximating the inverse of matrix T, it will now be shown how dual photography can be used with an immersive display system to achieve advanced and complex setups.

With reference to FIG. 11, in a preferred embodiment, a conventional front projector P1 (similar to projector 21 of FIG. 2) is used in conjunction with an immersive projector P2. As indicated by field-of-view lines 91, in the presently preferred embodiment the portion of a display surface covered by the field-of-view (i.e. FOV) of front projector P1 is a subset of the FOV of immersive projector P2, as indicated by field-of-view lines 93. Since the FOV of front projector P1 is entirely within the scope of the FOV of immersive projector P2, it would be desirable to have immersive projector P2 simulate the projected image produced by front projector P1. It should be emphasized, however, that in general, the FOV of front projector P1 does not necessarily need to overlap any part of the FOV of immersive projector.

Although, it is not necessary that the FOV of immersive projector P2 overlap any part of the FOV of front projector P1, it is desired that two transport matrices separately associating a camera C to front projector P1 and to immersive projector P2 be created. As it would be understood, the two transport matrices may be generated separately since the FOV's of P1 and P2 do not necessarily overlap.

However, in the specific example of the presently preferred embodiment, camera C is placed such that the FOV of camera C is a superset of the FOV of front projector P1 and a subset of the FOV of immersive projector P2. As indicated by field-of-view lines 95, the FOV of camera C completely encompasses the FOV of front projector P1, but is entirely engrossed by the FOV of immersive projector P2. To simulate the projected image from front projector P1 using immersive projector P2, one first determines a first light transport matrix, T₁, relating a first projected image p₁ from front projector P1 to a first captured image c₁ captured by camera C such that c₁=T₁p₁, and determines a second light transport matrices, T₂, relating a second projected image p₂ from immersive projector P2 to a second captured image c₂ captured by camera C such that c₂=T₂p₂. Consequently, one has the following relation c ₁ =T ₁ p ₁ and c ₂ =T ₂ p ₂ In order to simulate projected image p₁ from front projector P1 using immersive projector P2, one needs c₁ (i.e. the captured, projected image from front projector P1) to be the same as c₂ (i.e. the captured, projected image from immersive projector P2), i.e. one needs c ₂ =c ₁ which lead to the relation: T ₂ p ₂ =T ₁ p ₁ solving for p₂, one obtains the following relation: p ₂=(T ₂ ⁻¹)(T ₁ p ₁)

This means that to create image p₁, one can use front projector P1, or the same effect can be achieved by projecting the transformed image (T₂ ⁻¹)(T₁p₁) on immersive projector P2. This is shown in FIG. 12, where immersive projector P2 is used to simulate a front projector, such as projector P1 of FIG. 11. In FIG. 12, a virtual projector P1″, as simulated by immersive projector P2, is denoted by dotted lines. Therefore, image p₁, as projected by front projector P1 in FIG. 11, can be recreated by projecting transformed image (T₂ ⁻¹)×(T₁p₁) on immersive projector P2 in FIG. 12 to simulate a front projected image such that viewers 100 a, 100 b, and 100 c do not have to concern themselves with occluding any front projector.

An example of an image generated using this virtual projector implementation is shown in FIG. 13. A front projector image 101 is simulated using a large field-of-view display system. Projector 103, located along the bottom of FIG. 13 is part of the large field-of-view display system, and is used to generate image 101 shown in the center of FIG. 13.

FIGS. 14A to 14C illustrate the quality of the simulation by showing the real and simulated displayed images side-by-side. FIG. 14A shows the right side of an image projected by a real front projector, such as P1 of FIG. 11. FIG. 14B shows the corresponding left side of the image shown in FIG. 14A, but in FIG. 14B the left side of the image is projected by an immersive projector, such as P2 of FIG. 12, to simulate a virtual front projector, such as P1″ of FIG. 12. The quality of the simulated front projection, as created by the immersive projector is better illustrated in FIG. 14C, where the right side image of FIG. 14A is shown joined to the left side image of FIG. 14B, side-by-side, resulting in a seamless registration of the images created by a real front projector and a simulated, virtual front projector.

Two additional examples showing side-by-side comparisons of images created by a real front projector and images created by an immersive projector to appear as if they were generated by a front projector are shown in FIGS. 15A and 15B. In both FIGS. 15A and 15B, the left half of the shown image is created by an immersive projector to simulate a display from a virtual front projector, and the right side half of the shown image is created by a real front projector.

An alternate application of the present technique is better understood with reference to FIGS. 16A to 16B. In the present example, immersive projector P2 of FIG. 16C will be used to create various ambient lighting effects. To achieve this, camera C is positioned such that its FOV covers a significant portion of a display room 111, as shown in FIG. 16A. In FIG. 16A, camera C and immersive projector P2 are positioned such that the FOV of camera C encompasses most, if not all of the FOV of immersive projector P2. In the present example, P2 is shown as an immersive projector, but projector P2 may be any type of a projector, such a front projector. To establish a relationship between camera C and projector P2, a light transport matrix T₃ relating camera C to projector P2 is captured, i.e. determined, using any of the methods described above. Once this is done, an image c₃ as viewed (i.e. captured) by camera C will be related to a projected image p₃, as projected by projector P2, according to the following relationship: c ₃ =T ₃ p ₃ which results in p ₃=(T ₃ ⁻¹)×(c ₃) Consequently, one can build a virtual model of display surfaces of room 111. The virtually modeled room 111″, shown in FIG. 16B, may be a computer simulation, for example. Once virtual room 111″ is created, various simulated lighting effects may be added to modeled room 111″. For example, FIG. 16B shows virtual model room 111″ being lit by candle light from a large candle 113. The computer model further models the position and resolution of camera C (of FIG. 16A), shown as dotted box C in FIG. 16B. The computer model then “captures” (i.e. creates) a synthetic view c₃″ of simulated room 111″ from the viewpoint camera C to simulate a real image of virtual room 111″ as if it had been taken by real camera C of FIG. 16A. The simulated lighting effects of FIG. 16B can then be recreated in real room 111 of FIG. 16C using P2 by projecting the transformed image (T₃ ⁻¹)×(c₃″).

An example of an application of this technique is shown in FIG. 17. In the present case, it is desired to project an image 117 that is bigger than the wall of a real room 111. As was discussed above, various techniques may be used to calibrate a real projector to compensate for the angles of the walls and ceiling to the projection wall of real room 111, but the present invention solves this problem using a different approach. In the present example, virtual room 111″ (of FIG. 16B) has dimensions similar to real room 111 (of FIGS. 16C and 17), and image 117 is superimposed (i.e. projected) onto virtual room 111″. Image distortions, if any, can be readily corrected within the modeled room 111″. A synthetic view c₃″ of image 117 without distortion on simulated room 111″ from the viewpoint camera C is then created. Immersive projector P2 is then made to project transformed image (T₃ ⁻¹)×(c₃″) to recreate the corrected oversized image 117 on a wall of real room 111. As is shown in FIG. 17, the result is an undistorted projection that did not require calibrating projector P2 to compensate for curvatures (or other irregularities) on a projection surface.

The above described techniques may be applied to the creation of large field-of-view, i.e. FOV, displays. A larger FOV increase the sense of immersion and provides a more engaging experience for a viewer. The present approach describes an immersive projection system with a very large FOV. The system is also designed with a built-in large FOV camera/light sensor that is able to capture light from the areas covered by the projection FOV. The sensor allows the system to adapt the projected light so as to optimize image quality and more generally allow the system to interact with its environment. Although the present system is primarily motivated by the desire to display surround video content, it is important to note that this new projection system can also be used to view conventional video content.

With reference to FIG. 18, an exemplary projection system in accord with the present invention in its minimal form consists of the following components: a projector 121; a camera 123, which can be a digital still camera or a video camera; curved mirror 125, which can be spherical or otherwise; and mounting mechanisms for the above components. Light from projector 121 is reflected off curved mirror 125 before reaching a display surface 127, which can be any surface, including building walls, floors, ceilings, and dedicated projection screens. Display surface 128 can also be arbitrarily shaped. Reflecting the projected light off the curved mirror enlarges the projector FOV. Light rays from the environment, which may or may not have originated from the projector, also reflect off the curved mirror 125 before reaching the camera. This similarly enlarges the camera FOV.

FIG. 19 shows a prototype based on the design of FIG. 18, and all elements in FIG. 19 similar to those of FIG. 18 have similar reference characters and are described above. In this prototype, curved mirror 125 is hemispherical, and camera 123 is a high resolution digital still camera. FIG. 20 is an alternate view of the setup of FIG. 19, shows the view of mirror 123 as seen (very roughly) from camera 123. As can be seen in FIG. 20, camera 123 is able to see the floor, at least three vertical walls, and the ceiling from the reflection in mirror 125.

In FIG. 21A a room with the present projection system installed is shown under ambient lighting, and the same room is shown in FIG. 21B under immersive projection lighting from the present projection system. As can be seen, the present projection system is able to project images onto the two walls as well as the ceiling.

In the present embodiment, the FOV of projector 12 and the FOV of camera 123 are in general different, and may or may not overlap. When there is a significant overlap in the two FOV's, images captured by camera 123 can be used as feedback for improving the quality of a projected image from projector 121 in a manner similar to those described above. For example, feedback from to projector 121 can be used to compensate for variations in the display surface reflectance properties and shape so that the projected image appears as though it was projected on a flat white surface.

The FOV of camera 123 camera may also include areas not covered by the FOV of projector 121. For example, while projector 121 covers the front and side walls of the test room shown in FIGS. 21A and 21B, the camera may capture areas outside the projector FOV, possibly including areas where viewers are located. This allows the system to adapt and interact with viewers either by detecting and tracking the viewers or pointing devices. It may be possible for camera 123 to track small lights mounted, for example, on remote controls and facilitate user interaction.

In several of the designs described above, the projector and cameras do not have common optical centers. However it is possible to design projector-camera pairs with collocated optical centers. While a system with collocated optical centers allows the projector and the camera to have perfectly identical FOV, a system with non-collocated optical centers has the potential to allow 3D reconstruction of its environment.

With reference to FIG. 22, an alternative design for ceiling-mounted operation is shown. All elements similar to those of FIG. 18 have similar reference characters and are described above.

In FIGS. 23 and 24, two additional alternate configurations are shown. All elements in FIGS. 23 and 24 similar to those of FIG. 18 have similar reference characters and are described above. In FIG. 23 a planar mirror 141 is used to fold the optical path so that projector 121 and camera 123 can be placed under the curved mirror 125, thereby achieving a smaller footprint. FIG. 24 shows a booth design for enclosing projector 121, camera 123, curved mirror 125, and flat mirror 141 within a booth 143 for display booth operation.

Two additional variations on the present design are further shown in FIGS. 25A and 25B. FIG. 25A shows a design that uses a single curved mirror 125 and multiple projector-camera pairs 145. FIG. 25B shows a design that uses a single mirror pyramid 151 and multiple projector-camera pairs 145. With the latter, the optical centers of all projectors can be collocated within the mirror pyramid, creating a single virtual large FOV projector. Similarly, the camera optical centers can also be collocated to create a single virtual large FOV camera.

Finally, FIG. 26 shows that multiple large FOV projectors 153 a and 153 b can be used to achieve an even larger overall projection FOV. One or more conventional projectors 155 can also be used in combination. As seen, the FOV of projector 153 a, as indicated by dotted line 157 overlaps the FOV of projector 154 b, as indicated by dotted line 159, by an overlap amount 161.

Color-Based Feature Identification and Image System Calibration

Described above were systems and methods for obtaining transformations between projector and camera coordinates. One of the key measures of performance in calibrating a projector-camera system is the amount of time needed for capturing a light transport matrix, which matrix may be used to calibrate the system. Typically, the amount of time needed for capturing a light transport matrix is typically correlated with the number of images that need to be captured. Disclosed herein are systems and methods for using color to help identify feature points, which can be used to help reduce the number of display images used in calibrating an image system, such as a projector-camera system.

In the following description, for purpose of explanation, specific details are set forth in order to provide an understanding of the invention. It will be apparent, however, to one skilled in the art that the invention may be practiced without these details. One skilled in the art will recognize that embodiments of the present invention, some of which are described below, may be incorporated into a number of different systems and devices. The embodiments of the present invention may be implemented in software, hardware, firmware, or combinations thereof. Components shown below in block diagrams are illustrative of exemplary embodiments of the invention and are meant to avoid obscuring the invention. It shall also be understood that throughout this discussion that components may be described as separate functional units, which may comprise sub-units, but those skilled in the art will recognize that the various components, or portions thereof, may be divided into separate components or may be integrated together, including integrating within a single system or component. Furthermore, the components may comprise additional functionality, including additional functionality to support the disclosed functionality.

Connections between components within the figures are not intended to be limited to direct connections. Rather, data between these components may be modified, re-formatted, or otherwise changed by intermediary components. Also, additional or fewer connections may be used. It shall also be noted that the terms “coupled,” “communicatively coupled,” “receives as an input,” “coupled to receives information,” or the like shall be understood to include direct connections, indirect connections through one or more intermediary devices, and wireless connections.

Reference in the specification to “one embodiment” or “an embodiment” means that a particular feature, structure, characteristic, or function described in connection with the embodiment is included in at least one embodiment of the invention but may be in more than one embodiment. Also, the appearances of the phrase “in one embodiment” or “in an embodiment” in various places in the specification are not necessarily all referring to the same embodiment. It shall be noted that a set or a group need not contain a plurality of members of the set or the group. It shall also be noted that a point may comprise one or more pixels.

FIG. 27 depicts an example of an image system 2700 according to embodiments of the invention. Illustrated in FIG. 27 is a projector-camera system 2700, which comprises a camera 2710 and a projector 2715, each of which is communicatively coupled to a computing system 2705. It should be noted that although computing system 2705, camera 2710, and projector 2715 of the projector-camera system 2700 are depicted as separate components, one or more of those components may be integrated into a single device. It should also be noted that an image system may contain one or more cameras and/or one or more projectors. Other examples of image systems that may be used include, without limitation, the projector systems previously depicted or discussed. It should also be noted that the present invention may be used to calibrate one or more cameras or one or more projectors, in which one or more projectors or cameras may be removed or added to the image system. In embodiments, additional calibration may be performed for added devices.

A. Computing System Implementations

Aspects of the present invention may be implemented in any device or system or combination of devices or systems capable of performing according to the teachings of the present invention. Such devices and/or systems include, but are not limited to, a computer or a device with one or more processors, including without limitation, a general-purpose computer and a specific computer, such as one intended for calibrating image systems. The present invention may also be implemented into other computing devices, including without limitation, a digital camera, a projector, a multimedia device, and any other device that is capable of processing data. Furthermore, within any of the devices, aspects of the present invention may be implemented in a wide variety of ways including software, hardware, firmware, or combinations thereof. For example, the functions to practice various aspects of the present invention may be performed by components that are implemented in a wide variety of ways including discrete logic components, one or more application specific integrated circuits (ASICs), and/or program-controlled processors. It shall be noted that the manner in which these items are implemented is not critical to the present invention.

FIG. 28A depicts a functional block diagram of a system 2850 that may be used for feature identification using color information according to embodiments the present invention. As depicted, system 2850 comprises a calibration processor 2854 communicatively coupled to an image generator 2852 and a captured image processor 2856. In embodiments, system 2850 is computing system 2705.

In embodiments, image generator 2852 causes a calibration image to be displayed. In embodiments, the calibration image may be formed from a plurality of display images. The display image or images are configured to form a calibration image comprising a set of features and color patterns. Each feature in the set of features may be associated with a color pattern that identifies the feature. In embodiments, a feature may be formed from feature components; for example, two lines may be feature components and their intersection may form a feature. Because features can be uniquely identified using color information in the calibration image, the number of display images that would typically be required to form a calibration image may be reduced by displaying a plurality of feature or feature component in a display image. That is, because color information can be used to identify features, multiple feature or feature component may be display in a single display image.

Captured image processor 2856 receives at least a portion of the calibration image. In embodiments, captured image processor 2856 may receive from a camera captured images of at least portions of the display image or images. Because the field of view of the capturing device, whether integrated with captured image processor 2856 or separate but interfaced with captured image processor 2856, may not cover the entire display image, only a portion of a display image may be captured. If only a portion of the display image or images is captured, only a portion of the calibration image may be formed. Captured image processor 2856 forms the calibration image from the captured display image or images. In embodiments, captured image processor 2856 may form the calibration image by overlaying the captured display images. In embodiments, captured image processor 2856 identifies at least some of the features in the at least a portion of the calibration image using color patterns in the at least a portion of the calibration image. Examples of how color patterns may be used to identify features is described in more detail below.

Information about the identified features, such as the position information of the features in the captured image, is provided to a calibration processor 2854. Because calibration processor 2854 also receives information from image generator 2852 regarding the positions of the features in the displayed image, calibration processor 2854 may use that information to calibrate between the projector used to display the display images and the camera used to capture the display images. Calibration processor 2854 may use homography transform techniques as discussed above in calibrating to obtain a light transport matrix. In embodiments, calibration processor 2854 may also calculate an inverse light transport matrix for use in an image system to calibrate multiple projectors, as described previously.

In embodiments, the component of system 2850 may be executed on a instruction-executing/computing device. FIG. 28B depicts a functional block diagram of an embodiment of an instruction-execution/computing device 2875 that may implement or embody embodiments of the present invention. As illustrated in FIG. 28B, a processor 2802 executes software instructions and interacts with other system components. In an embodiment, processor 2802 may be a general purpose processor such as (by way of example and not limitation) an AMD processor, an INTEL processor, a SUN Microsystems processor, or a POWERPC compatible-CPU, or the processor may be an application specific processor or processors. A storage device 2804, coupled to processor 2802, provides long-term storage of data and software programs. Storage device 2804 may be a hard disk drive and/or another device capable of storing data, such as a computer-readable media (e.g., diskettes, tapes, compact disk, DVD, and the like) drive or a solid-state memory device. Storage device 2804 may hold programs, instructions, and/or data for use with processor 2802. In an embodiment, programs or instructions stored on or loaded from storage device 2804 may be loaded into memory 2806 and executed by processor 2802. In an embodiment, storage device 2804 holds programs or instructions for implementing an operating system on processor 2802. In one embodiment, possible operating systems include, but are not limited to, UNIX, AIX, LINUX, Microsoft Windows, and the Apple MAC OS. In embodiments, the operating system executes on, and controls the operation of, the computing system 2875.

An addressable memory 2806, coupled to processor 2802, may be used to store data and software instructions to be executed by processor 2802. Memory 2806 may be, for example, firmware, read only memory (ROM), flash memory, non-volatile random access memory (NVRAM), random access memory (RAM), or any combination thereof. In one embodiment, memory 2806 stores a number of software objects, otherwise known as services, utilities, components, or modules. One skilled in the art will also recognize that storage 2804 and memory 2806 may be the same items and function in both capacities.

In an embodiment, computing system 2875 provides the ability to communicate with other devices, other networks, or both. Computing system 2875 may include one or more interfaces or adapters 2812, 2814 to communicatively couple computing system 2875 to other networks and devices, such as camera 2710 and projector 2715. For example, computing system 2875 may include a network interface 2812, a communications port 2814, or both, each of which are communicatively coupled to processor 2802, and which may be used to couple computing system 2875 to other computer systems, networks, and devices.

In an embodiment, computing system 2875 may include one or more output devices 2808, coupled to processor 2802, to facilitate displaying graphics and text. Output devices 2808 may include, but are not limited to, a display, LCD screen, CRT monitor, projector, printer, touch screen, or other device for displaying information. Computing system 2875 may also include a graphics adapter (not shown) to assist in displaying information or images on output device 2808.

One or more input devices 2810, coupled to processor 2802, may be used to facilitate user input. Input device 2810 may include, but are not limited to, a pointing device, such as a mouse, trackball, or touchpad, and may also include a keyboard or keypad to input data or instructions into computing system 2875.

In an embodiment, computing system 2875 may receive input, through communications port 2814, network interface 2812, stored data in memory 2804/2806, and/or through an input device 2810.

One skilled in the art will recognize no computing system is critical to the practice of the present invention. One skilled in the art will also recognize that a number of the elements described above may be physically and/or functionally separated into sub-modules or combined together.

It shall be noted that embodiments of the present invention may further relate to computer products with a computer-readable medium that have computer code thereon for performing various computer-implemented operations. The media and computer code may be those specially designed and constructed for the purposes of the present invention, or they may be of the kind known or available to those having skill in the relevant arts. Examples of computer-readable media include, but are not limited to: magnetic media such as hard disks, floppy disks, and magnetic tape; optical media such as CD-ROMs and holographic devices; magneto-optical media; and hardware devices that are specially configured to store or to store and execute program code, such as application specific integrated circuits (ASICs), programmable logic devices (PLDs), flash memory devices, and ROM and RAM devices. Examples of computer code include machine code, such as produced by a compiler, and files containing higher level code that are executed by a computer using an interpreter.

B. Color-Based Calibration Methods

Systems and methods for developing a homography transformation between projector and camera coordinates were previously described. For example, in embodiments, two sets of stripes at known projector x (vertical stripes) and y (horizontal stripes) coordinates may displayed. In order to identify each stripe uniquely, each stripe is sequentially projected (i.e., one stripe after another). Thus, at any time, at most one stripe is projected. If more than one stripe may be projected and identified at a time, it would reduce the number of images that need to be captured and thereby reduces the amount of time needed for the calibration process.

FIG. 29 illustrates a color-based method for calibrating an image system according to embodiments of the invention. In embodiments, a projector generates an image comprising a set of features and color patterns, each feature in the set of features being associated with a color pattern that identifies the feature. The image may be generated by displaying (2905) features with color information. As illustrated in the embodiments presented below, the features and color information of the calibration image may be obtained from more than one display image or may be obtained from a single display image.

The displayed information is captured (2910) by one or more cameras. The captured image or images contains color information that is used to uniquely identify (2915) at least some of the features. Having identified at least some of the features, the feature information may be used to calibrate (2920) an image system. In embodiments, the image system, such as system 2700, may be calibrated using one or more of the methods discussed previously. For example, a homography transform may be used to obtain a light transport matrix.

In embodiments, features may be obtained by projecting simultaneously features or feature components rendered with distinct colors. One possible color palette is shown in FIG. 30. FIG. 30 depicts a set of colors that may be used in a color-based method for calibrating an image system according to embodiments of the invention. This particular palette comprises the corners of the Red-Green-Blue color cube, with the exception of the black corner (0, 0, 0).

For identifying the respective colors of each projected stripe, color matching systems and methods, which are known to those skilled in the art, may be used. Additional examples of color matching includes, by way of illustration and not limitation, the systems and methods disclosed in commonly-owned: U.S. Pat. No. 6,038,340, entitled “System and method for detecting the black and white points of a color image,” filed on Nov. 8, 1996, listing Hakan Ancin, Anoop Bhattacharjya, and Joseph Shu as inventors; U.S. Pat. No. 6,185,004, entitled “Self-calibration for color image reproduction system,” filed on Mar. 30, 1998, listing Tsung-Nan Lin and Anoop Bhattacharjya as inventors; and U.S. Pat. No. 6,373,590, entitled “Method and apparatus for slant adjustment and photo layout,” filed on Feb. 4, 1999, listing Hakan Ancin, Sean Miceli, Sean Gugler, and Anoop Bhattacharjya as inventors. Each of the above-listed patents is incorporated herein by reference in its entirety.

1. Color Lines

FIG. 31 illustrates a color lines method for calibrating an image system according to embodiments of the invention. In embodiments, using a color palette such as the one depicted in FIG. 30, seven lines may be simultaneously projected (3105) because each of the lines may be identified by its unique color.

For example, a first set of colored lines or stripes may be projected (3105). FIG. 32A illustrates a display of horizontal color lines that may be used for calibrating an image system, such as a projector-camera system, according to embodiments of the invention. The first set of uniquely-colored color lines are captured (3110) by a camera. A second set of uniquely-colored lines that are non-parallel to the first set of lines may then be projected (3115). FIG. 32B illustrates a display of vertical color lines that may be used for calibrating an image system according to embodiments of the invention. After the second set of lines are captured (3120) by the camera, feature information may be obtained. The intersections of the horizontal and vertical feature components create features that are identifiable. For example, there is only one feature that is formed by the intersection of a horizontal red line and a vertical green line. By projecting two images of seven stripes, up to 49 feature points may be identified. By way of comparison, using monochromatic stripes, 14 images would have to be captured to identify the same number of features. One skilled in the art shall recognize that more features may be captured by having more entries in the color palette and thereby displaying more x-lines and y-lines. Because the feature locations in the projector coordinates are known and the corresponding features in the camera coordinates may be identified using the color information, a homography transform may be used to calibrate (3125) the system.

2. Color Grid

In an alternative embodiment, a color grid may be used for capturing projection geometry. FIG. 33 illustrates a color grid method for calibrating an image system according to embodiments of the invention. Instead of projecting vertical and horizontal lines and using their intersection to identify feature points, a color grid may be projected (3305) and captured (3310). The color coding of the grid may be unique so that each feature point is distinguishable from another. FIG. 34A illustrates an example of a color grid using colors from the color palette depicted in FIG. 30.

In FIG. 34A, the top row and leftmost column of the table shows the palette colors, and the table shows how the colors are assigned to each grid. In the grid, each grid corner can be uniquely identified by the color combination in the vicinity of the grid. For example, feature 3405 may be uniquely identified as the point that has Red at the upper left of the feature point, Yellow at the upper right of the feature, Yellow at the lower left of the point, and Red at the lower right of the feature. No other feature point in the grid has that same color pattern configure. Consider, for example, feature 3410 that has Green at the upper left, White at the lower left, and black at the right.

It should be noted that the color palette can be different, and the order of the colors can also be permuted to produce a different color grid design. For example, the design in FIG. 34B is different from the grid shown in FIG. 34A, but they both have unique color combinations at each grid point. Thus, one skilled in the art shall recognize that no particular color palette or color grid pattern is critical to the present invention.

Having displayed and captured the color grid, it is now possible to compute the homography transform (3315) by directly fitting a color grid template to the captured grid.

3. Line and Color Grid

In embodiments, a combination of monochromatic lines may be used to create feature points and a color grid may be used to identify each point uniquely. FIG. 35 illustrates an alternative color grid method 3500 for calibrating an image system according to embodiments of the invention.

As depicted in FIG. 35, a first set of lines may be displayed by a projector and captured (3505) by a camera. FIG. 36A illustrates an example of horizontal lines, which may be monochromatic lines that are used in conjunction with a color grid for calibrating a projector-camera system according to embodiments of the invention. A second set of lines may next be displayed by a projector and captured (3510) by a camera. FIG. 36B illustrates an example of vertical lines, which may be monochromatic lines that are used in conjunction with a color grid for calibrating a projector-camera system according to embodiments of the invention.

FIG. 36C illustrates features obtained from the intersection of the feature components displayed in FIGS. 36A and 36B. These features may be used in conjunction with a color grid for calibrating an image system according to embodiments of the invention. The intersections of the horizontal lines (3610-1 through 3610-6) and the vertical lines (3615-1 through 3615-8) create a set of 48 features (3620-1 through 3620-48) as show in FIG. 36C. One skilled in the art shall recognize that a different number of lines and features or a different configuration of features may be used.

In this embodiment, because the horizontal and vertical feature components were monochromatic, the features (3620) may be unambiguously identifying by projecting (3515) a color grid. FIG. 36D illustrates an overlay of the features obtained from the intersection of the lines displayed in FIGS. 36A and 36B and a color grid that may be used for calibrating an image system according to embodiments of the invention.

To assign projector coordinates to each of the 48 points, the location of each of the points are overlaid on the projected color grid, as shown in FIG. 36D. Since the lines projected correspond exactly to the edges on the color grid, the feature points detected correspond to the corners on the grid. Using the color combination in the vicinity of each feature, it is possible to assign (3520) the correct projector coordinates to each feature 3620 since the corresponding projector coordinates on the grid are known. Using some of all of the feature information, the projector system may be calibrated (3525).

While the invention has been described in conjunction with several specific embodiments, it is evident to those skilled in the art that many further alternatives, modifications, and variations will be apparent in light of the foregoing description. Thus, the invention described herein is intended to embrace all such alternatives, modifications, applications and variations as may fall within the spirit and scope of the appended claims. 

1. A method for identifying predefined point locations within a captured-image, said method comprising the following step: (a) defining a calibration image having a plurality of said predefined point locations distributed within it, each of said plurality of predefined point locations being uniquely identified and distinguished from each other within the calibration image by a unique combination of adjacent colors around each predefined point location; (b) generating said captured-image by capturing at least a portion of said calibration image; and (c) identifying point locations within the captured-image that correspond to predefined point locations in said calibration image according to their matching unique combination of adjacent colors, correlating the identified point locations within the captured-image to their corresponding predefined point locations within said calibration image; wherein in step (a), said calibration image includes: a first group of colored lines arranged in parallel to each other, each colored line being distinct from all other lines in said first group by color; a second group of colored lines arranged in parallel to each other, each colored line being distinct from all others lines in said second group by color, said second group of colored lines being arranged to intersect said first group of colored lines, wherein said predefined point locations are defined by the intersections of color lines from said first group with color lines from said second group, each of said intersections being defined by a unique combination of intersecting colors.
 2. The method of claim 1, wherein in step (b), the capturing of at least a portion of the calibration image is by means of a digital camera.
 3. The method of claim 2, further comprising: after step (a) and before step (b), projecting said calibration image by means of a digital projector, and correlating each of said predefined point locations within said calibration image with their corresponding array locations within the projector pixel array of said digital projector; and wherein in step (b), said digital camera captures the projections of said digital projector.
 4. The method of claim 3, wherein said digital projector projects said calibration image in its entirety as a single projected image, and said digital camera captures said at least a portion of said calibration image by capturing at least a portion of the area of said calibration image as projected in its entirety in said single projected image.
 5. The method of claim 3, further including after step (c): correlating the identified point locations within the captured-image to array locations within the digital projector's projector pixel array according to their commonly correlated predefined point locations within said calibration image.
 6. The method of claim 3, further including after step (c): applying a homography transform to calibrate said digital projector to said digital camera according to their mutually corresponding predefined point locations within said calibration image.
 7. The method of claim 1, wherein the number of predefined point locations is limited by the number of available unique color combinations in said calibration image.
 8. The method of claim 1, wherein said calibration image has more than 28 of said predefined point locations individually distinguishable from each other.
 9. The method of claim 1, wherein each colored line in said first and second groups is a solid color.
 10. The method of claim 1, wherein each predefined position is defined by a set of four quadrant areas joined at each respective predefined position, each set of four quadrant areas conforming to the four quadrants of a Cartesian plane with their adjoining predefined position at the Cartesian origin, each quadrant area within each set of four quadrant areas having a distinct color and each set of four quadrant areas forming a unique color pattern distinct from all other sets of four quadrant areas.
 11. The method of claim 10, wherein each of said quadrant areas has a square shape and the arrangement of said predefined point positions forms a checkered pattern on said calibration image, with each cluster of four squares adjoining a predefined position forming a color pattern distinct from all other clusters of four squares adjoining all other predefined positions within said calibration image.
 12. A non-transitory computer readable medium comprising instructions for implementing the method of claim
 1. 13. A method of calibrating a digital camera to a digital projector, comprising: (a) defining a calibration image having a plurality of predefined point locations distributed within it, each of said plurality of predefined point locations being uniquely identified and distinguished from each other within the calibration image by a unique combination of adjacent colors around each predefined point location; (b) having said digital projector project a projected-image representation of said calibration image; (c) having said digital camera capture the projections of said digital projector; (d) defining a captured-image representation of said calibration image from the output of said digital camera; and (e) identifying captured point locations within the captured-image representation of said calibration image that correspond to predefined point locations within the projected-image representation of said calibration image according to their matching unique combination of adjacent colors, and correlating the identified captured point locations within the captured-image representation of said calibration image to their corresponding predefined point locations within the projected-image representation of said calibration image; wherein in step (a), said calibration image includes: a first group of colored lines arranged in parallel to each other, each colored line being distinct from all other lines in said first group by color; a second group of colored lines arranged in parallel to each other, each colored line being distinct from all others lines in said second group by color, said second group of colored lines being arranged to intersect said first group of colored lines, wherein said predefined point locations are defined by the intersections of color lines from said first group with color lines from said second group, each of said intersections being defined by a unique combination of intersecting colors.
 14. The method of claim 13, wherein: in step (b), said digital projector creates said projected-image representation of said calibration image by creating a projection of said calibration image in its entirety as a single projected image; and in step (c), said digital camera captures said single projected image as a single captured image; in step (d), said single captured image is defined as said captured-image representation of said calibration image.
 15. The method of claim 14, wherein each predefined position within said calibration image is defined by a set of four quadrant areas joined at each respective predefined position, each set of four quadrant areas conforming to the four quadrants of a Cartesian plane with their adjoining predefined position at the Cartesian origin, each quadrant area within each set of four quadrant areas having a distinct color and each set of four quadrant areas forming a unique color pattern distinct from all other sets of four quadrant areas within said calibration image.
 16. The method of claim 15, wherein each of said quadrant areas has a square shape and the arrangement of said predefined point positions forms a checkered pattern on said calibration image, with each cluster of four squares adjoining a predefined position forming a color pattern distinct from all other clusters of four squares adjoining all other predefined positions within said calibration image.
 17. A method of calibrating a digital camera to a digital projector, comprising: (a) defining a calibration image having a plurality of predefined point locations distributed within it, each of said plurality of predefined point locations being uniquely identified and distinguished from each other within the calibration image by a unique combination of adjacent colors around each predefined point location; (b) having said digital projector project a projected-image representation of said calibration image; (c) having said digital camera capture the projections of said digital projector; (d) defining a captured-image representation of said calibration image from the output of said digital camera; and (e) identifying captured point locations within the captured-image representation of said calibration image that correspond to predefined point locations within the projected-image representation of said calibration image, according to their matching unique combination of adjacent colors, and correlating the identified captured point locations within the captured-image representation of said calibration image to their corresponding predefined point locations within the projected-image representation of said calibration image; wherein: in step (b), said digital projector creates said projected-image representation of said calibration image by creating a projection of said calibration image in its entirety as a single projected image; in step (c), said digital camera captures said single projected image as a single captured image; in step (d), said single captured image is defined as said captured-image representation of said calibration image; each predefined position within said calibration image is defined by a set of four quadrant areas joined at each respective predefined position, each set of four quadrant areas conforming to the four quadrants of a Cartesian plane with their adjoining predefined position at the Cartesian origin, each quadrant area within each set of four quadrant areas having a distinct color and each set of four quadrant areas forming a unique color pattern distinct from all other sets of four quadrant areas within said calibration image; and prior to step (b), said digital projector creates a first projection of a first plurality of monochromatic lines followed by a second projection of a second plurality of monochromatic lines, said first plurality of monochromatic lines being arranged in parallel to each other, said second plurality of monochromatic lines being arranged in parallel to each other, all monochromic lines from said first and second pluralities of monochromatic lines being the same color, and the intersections of monochromatic lines from said first and second pluralities of monochromatic lines coincide with said predefined point locations.
 18. The method of claim 17 further including: extracting light transport coefficients from the intersections of the monochromatic lines of said first and second pluralities of monochromatic lines.
 19. A non-transitory computer readable medium comprising instructions for implementing the method of claim
 13. 